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ABSTRACT 



We report the discovery of strong soft X-ray emission lines and a hard continuum above 2 keV 
in the Narrow-Line Seyfert 1 galaxy Mrk 335 during an extremely low X-ray flux state. Mrk 335 
was observed for 22 ks by XMM-Newton in July 2007 as a Target of Opportunity to examine it 
in its X-ray low-flux state, which was discovered with Swift. Long-term light curves suggest that 
this is the lowest flux state this AGN has ever been seen in. However, Mrk 335 is still sufficiently 
bright that its X-ray properties can be studied in detail. The X-ray continuum spectrum is very 
complex and requires several components to model. Statistically, partial covering and blurred 
reflection models work well. We confirm the presence of a strong narrow Fe line at 6.4 keV. High- 
resolution spectroscopy with the XMM-Newton-KGS reveals strong, soft X-ray emission lines not 
detected in previous, higher signal-to-noise, XMM-Newton observations, such as: highly ionized 
Fe lines, O VII, Ne IX and Mg XI lines. The optical/UV fluxes are similar to those previously 
measured with Swift. Optical spectroscopy taken in 2007 September do not show any changes to 
optical spectra obtained 8 years earlier. 

Subject headings: galaxies: active, galaxies: individual (Mrk 335), galaxies: Seyferts, X-rays: galaxies, 
ultraviolet: galaxies 



1. Introduction 

The Narrow-Line Seyf ert I galaxy (NLSI; 
Osterbrock fe Poggej Il985l ) Mrk 335 (a 20 oo = 
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00 h 06 m I9. s 5, (S2000 = +20°I2'11'.'0, z=0.026) was 
discovered as a bright X -ray AGN with UHURU 
(jTananbaum et al.lll9781 ) in 1971. Consequently 
it has been a target of almost all X-ray observato- 
ries including: EINSTEIN, EXOSAT, GINGA, 



ern 



ROS AT, ASCA, XMM and Suzaku , 

19821 IPounds et al.lll987tlNandra fe Poundslll994l; 



Grup e et al.ll20 01: Geor ge et al.l l2000: Gondoi n et al 
20021: Longinotti et al. 2007a Longinotti et al 



l2007bl : lO'Neill et alJ 120071 : lLarsson et all l2008h . 

The X-ray spectrum of Mrk 335 appears to be 
rather complicated and cannot be modeled by 
a single power law spectrum. Possible interpre- 
tations for the X-ray sp ectrum include: absorp- 
tion by ionized gas (e.g. Nandra fe Pounds] 1994: 
| Leighlvlll999bh . partial-covering (e.g. lTanaka et al 
l2005j), an d X-ray reflection on the accretion disk 
"e.g. Ballantvne et al. 2 00 it I Crummy et al 



20061: lLonginotti et al.ll2007at lO'Neill et al ll2007t 
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Larsson et al.ll2008l ) . Support for the reflection in- 
terpretation comes from the pre sence of a broad- 



ened and asymmetric Fe Kg line (jLonginotti et al 
2007at lO'Neill et af] l2007h . Such lines are pre- 



dicted as a consequence of reflection of the pri- 
mary powe r law on the accretio n disk close to the 
black hole (Fabia n et al. and arc seen in a 

erowme; number of AGN (see iNandra et al. Il2007t 
and references therein). Mrk 335 is a highly vari- 
able X-ray source that typically var ies by factors 
of 2 to 4 on hourly time scales (e.g. 



19931 : lO'Neill et al.ll2007l) 



Turner et al 



However, when Mrk 335 was observed with 
the Gamma-Ray Bur st Explorer Mission Swift 
in 2007 May and June 
it appeared significantly 



Gehrels et al. 



Grupe et al 



2004) 
2007rJl 



fainter in its X-ray flux by factors of more than 
30 compared to the 2006 XMM-Newt o n and 
Suzaku observatio ns j Longinotti et al. 2007bt 



O'Neill et all l2007t lLarsson et alj 120081 ). This 



drop in X-rays coincided with substantial changes 
in the X-ray spectrum. Mrk 335 showed a very 
hard X-ray component at energies above 2 keV, 
very similar to what has been seen in other NLS1 
galaxi es (e.g. Gallo 2006) and specifically in N GC 



4051 (jGuainazzi et alJll998t bonti et al.ll2Q06h 



As discussed in Grupe et al. ( 2007bl ). one pos 



sible explanation for the sudden drop in the X-ray 
flux of Mrk 335 is a partial covering absorber that 
inters ects our line-of-sight w ith the X-ray emitting 



Grupe et al. (|2007bl) found that when cor- 



zone. 

recting for this absorber the intrinsic X-ray flux 
of Mrk 335 between the high- and low-flux states 
varied only by factors of 4-6. 

In order to investigate the nature of the 
low-state in more detail we initiated a Target- 
of-Opportunity (ToO) observation with XMM- 
Newton, which was executed on 2007 July 10. 



XMM-Newton ToO observation. 

The paper is organized as follows: in §2 we 
present the XMM-Newton and Swift observations 
and describe the data reduction, in §3 we present 
the results that are discussed in §4. Through- 
out the paper spectral indeces are denoted as en- 
ergy spectral indeces with F v cx v~ a . Luminosi- 
ties are calculated assuming a ACDM cosmology 
with Om=0.27, J7a=0.73 and a Hubble constant 
of Ho=75 km s _1 Mpc -1 corresponding to a lu- 
minosity distance D=105 Mpc. All errors are 90% 
confidence unless stated otherwise. 

2. Observations and data reduction 

2.1. XMM-Newton and Swift Observations 



component in the spectrum associated with the 
red-wing of a relativistic iron line. This wing is 
not required if the spectrum is modified by a par- 
tial covering absorber. The brightness of Mrk 335 
provides us with the rare opportunity to study an 
AGN in its low-state, but still with sufficient sig- 
nal to obtain well-exposed spectra in a relatively 
short observation. In this paper we report the 



XMM-Newton (jJansen et all l200ll ) observed 
Mrk 335 on 2007 July 10 for a total of 22 ks (Ob- 
sID 0510010701). A summary of the observations 
with each of the instruments on-board XMM- 
Newton is given in Table[TJ The E uropean Pho 



ton I maging Camera (EPIC) pn (jStriider et al 



20011 ) was operated in Large Window mode with 
the thin filter. This combination was chosen to 
avoid pileup in ca se the AGN re-brigh tened. The 



two EPIC MOS dTurner et alj l200lj ) were both 



operated in Full-Frame mode with the medium 
filters. High-resolution X-ray spectroscopy was 
performed using the two Reflection Grating Spec- 



trometers (RGS; I den Herder et al.l 120011 ) on- 
board XMM-Newton. Optical photometry was 
perfor med in 5 filters with the Optical Monitor 
(OM; iMason et al~l l200lh . The data are used 



Newton in 2000 and 2006 dGondoin et al. 


2002| 


Crummv et al.l2006t Longinotti et al.l2007a 


O'Neill et al. 


2007). Of interest with Mrk 335 is the broad 



to measure the optical-to- X-ray spectral energy 
distribution of Mkn 335 during the XMM-Newton 
observation. Due to slew problems at the be- 
ginning of the observations, V-filter observations 
were not obtained. All OM observations were 
performed in a science-user defined configuration 
with a 7 x 7 observing window. 

The XMM-Newton data were processed in the 
standard way using the XMMSAS version xmm- 
sas_20060628-1801-7.0.0. The EPIC pn data 
were checked for episodes of high particle back- 
ground. The background count rate was found 
low throughout the observation. The source X-ray 
photons in the EPIC pn were selected in a circular 
region with a radius of 1 . Likewise, background 
photons were selected from a nearby, source-free 
region with the same radius. Only single and dou- 



2 



ble events (PATTERN . le . 4) and single to quadru- 
ple events (PATTERN. le . 12) were selected for the 
pn and MOS data, respectively. The spectra were 
rebinned with grppha version 3.0.0 with 50 pho- 
tons per bin. The redistribution matrices and the 
auxiliary response files were created by the XMM- 
SAS tasks rmfgen and arfgen, respectively. RGS 
spectra and response matrixes were created by the 
standard RGS XMMSAS tool rgsproc. The RGS 
spectra were rebinned with 10 photons per bin us- 
ing grppha. Spectral fits to the EPIC pn and MOS, 
and RGS spectr a were perform ed with XSPEC 



(NAOC) on 2007 September 07. Two spectra 
of Mrk 335 were taken with 30 and 60 min ex- 
posures during clear weather conditions and 2 
seeing. The spectrograph was equipped with a 
600 lines per mm grating and the 2 slit resulting 
in a resolution of 5A in first order, covering the 
wavelength range between 4377 to 7120A. The op- 
tical data were reduced and analyzed in a standard 
way with the ESO Munich Image Data Analysis 
System (MIDAS, version 06FEBpll.O). All line 
measurements were performe d on the coadded 



spectr um after subtracting the iBoroson fc Green 



12.3. lx (|Arnaud l Il996). The OM data (Il992l) I Zw 1 Fe II template as described in 



were processed with the XMMSAS task omichain. 
The magnitudes and fluxes of Mrk 335 were taken 
from the source lists created by the omichain task. 
For the count rate to flux conversion we used the 
conversion factors given in the OM Calibration 
document XMM-SOC-CAL-TN-0019. 

The SwiftX-r&y Telescope (XRT; iBurrows et~al 



Grupe et al 



(|2004al ). For the H/3 line we sub- 
tracted the narrow component by constructing a 
narrow-line template from the [01111 5007 line as 



described in lGrupe et al. I (|2004ah . 



Results 



20051) and UV/Optical Telescope fUVOT: iRoming et al. I 3-1 " Tem P oral Behavior 



2005) 2007 May and June observat ions and the 
data r eductions are described in iGrupe et al. 
(|2007bl ). In order to compare the photometry in 



the OM with the Swift UVOT we selected 5 field 
stars with similar brightness in V as Mrk 335. The 
measured (uncorrected) magnitudes of these stars 
are listed in TableEl W e confirm the earlier result 
bv lGrupe et al. (|2007bl ) that in U the magnitudes 
in the OM and UVOT agree perfectly. However, 
the OM magnitude had to be adjusted by —0.10 
mag in B and +0.30 in the UVM2 filter. In addi- 
tion to the 2000 XMM-Newton observation, during 
the 2007 XMM-Newton ToO observation we also 
observed in UV Wl and UV W2 as listed in Ta- 
bled] In these filters the OM magnitudes have 
to be adjusted, compared to the Swift UVOT, by 
+0.78 mag and +0.15 mag for the UVW1 and 
UVW2 filters, respectively. The 5 stars also show 
that the OM and UVOT detectors are stable and 
no variability is seen in these 5 stars. This suggests 
that any variability detected in the Optical/UV 
in Mrk 335 is real. A summary of the UVOT cali- 
bratio n has recently been published bv lPoole et al.l 
(|2007h . 



3.1.1. Long-Term Light curves 

Figured] displays the observed 2.0-10.0 kcV, 0.2- 
2.0 keV, and UV (2256A) long-term light curves 
of Mrk 335 covering the time between 1971 and 
2007 July. The X-ray fluxes and references are 
listed in Tabled The UV fluxes were derived from 
the U V IUE and HST database by 
2006). To remain consistent with 



Dunn et al. 



Dunn et al. 



2006), the XMM-Newton OM and Swift UVOT 



2.2. Optical Spectroscopic Observation 

Optical spectra of Mrk 335 were obtained with 
the 2.16 m telescope at the Xinglong Observatory, 
National Astronomical Observatories of China 



light curves were created in the UV M2 filters, 
since its cent ral wavelength is clo sest to 2256A. 
As noticed bv lPounds et all l| 19871 ), the EXOSAT 
data of Mrk 335 suggest that it was in an un- 
usual X-ray low-state in 1983. This seems coin- 
cident with a minimum in the UV light curve as 
well. However, about a year later in December 
1984, Mrk 335 was back in a typical high state. 
The plots in Figure[Tj demonstrate that from 2007 
May to July Mrk 335 was seen in an historic X- 
ray low-state. This drop happened in less than 
one year since Mrk 335 was last observ ed in June 
2006 with Suzaku ([Larsson et al.l 120081 ). In order 
to determine when and how fast this recent flux 
drop occurred, the RXTE ASM light curve was 
examined. The month-to-month ASM light curve 
indicates that Mrk 335 has been variable in X-rays 
over the last decade. However, Mrk 335 is too faint 
to constrain the occurrence of the most recent flux 
drop with the ASM. Mrk 335 also appears to be 
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fainter in the UV in 2007 compared with previ- 
ous epochs. Tableg] lists the XMM-Newton OM 
and Swift UVOT measurements. During the 2007 
May 17 Swift observation Mrk 335 appears to be 
fainter in the UV by about 0.2 mag compared to 
the Swift and XMM-Newton observations before 
and after. 

3.1.2. Short-term Variability 

The left panel of Figure[2] shows the 2007 July 
XMM-Newton EPIC pn light curve in the 0.2-10.0 
keV band. During the 16.8 ks exposure the AGN 
showed a steady decay in flux by about 15%. The 
light curve is relatively quiescent compared to the 
high-flux state variability seen in th e 2006 XMM- 



Newt o n and Suzaku light curves (|Q'Neill et al 
l2007t lLarsson et all I2008L respectively). The 
variability is in agreement with the variability - 
flux d ependence reported by lUttlev fe McHardv I 
(|200ll ). However, the 2007 July XMM-Newton 



observation was rather short (17 ks). During the 
low-state Swift observation on 2007 May 25, which 
covers a time period of 70 ks, Mrk 335 is quite vari- 
able by a factor of about 2 within 20 ks (right panel 
in Figure[2]). In both cases the AGN is variable 
and deviates significantly from a constant value. 
X 2 tests show that the \ 2 l v is 49/15 and 55/11 for 
the XMM-Newton pn and Swift XRT light curves, 
respectively, assuming the mean count rates as 
constant values. 

In Figure[3] the XRT and UVW2 light curves 
for all Swift segments between 2007 May 17 and 
July 02 are displayed (the magnitudes of segments 
001-005 are listed in Tabled]). The light curves 
show that Mrk 335 is quite variable on interme- 
diate time scales of days to weeks. Note that the 
UVW2 light curve follows the XRT count rate sug- 
gesting some connection between the Optical/UV 
and X-ray variability. The hardness ratio3 vari- 
ability curve may suggest some spectral variabil- 
ity in the X-ray band (middle panel of Figure[3]), 
however a x 2 test shows that this is not significant 
considering the uncertainties. 



1 The hardness ratio is denned as HR = (H - S)/(H + S) 
(where S and H are the number of photons in the 0.3-1.0 
keV and 1.0-10.0 keV bands, respectively). 



3.2. X-ray Spectral analysis 

Table[5] summarizes the results from the X-ray 
spectral analysis in the 0.3-12.0 keV range of the 
2007 XMM-Newton ToO (low-flux state) and 2006 
XMM-Newton GO (high-flux state) observations. 
The 2007 EPIC pn and MOS spectra cannot be 
fitted by a single absorbed power law model. This 
is similar to t he previous XMM-N ewton and Swift 
observations (jGrupe et al. Il2007bl) which cannot be 
fitted by a single power law model either. Ap- 
plying a broken power law model to the pn and 
MOS data significantly improves the fit, but as 
shown in Figure|4]the X-ray spectrum during the 
2007 XMM-Newton ToO observation is quite com- 
plicated and requires further components. There 
are strong residuals at energies below 1 keV and 
between 3-7 keV. The residuals at lower energies 
can be associated with strong X-ray emission lines 
found in the RGS spectra of Mrk 335 (see §EE23). 
Figure|3]also shows that the pn and MOS data are 
consistent with each other. Due to some internal 
processing problems, MOS data were not made 
available to us in time to be handled in the full 
analysis. In order to represent soft X-ray emission 
lines we added two Gaussian profiles at 0.5 and 
0.9 keV. These Gaussians do not represent spe- 
cific lines, but are used to mimic line emission and 
improve the fit. The emission feature at 6.4 keV is 
associated with near-neutral Fe Ka emission (e.g 
Fe I- XVII), and is modeled with a single, narrow 
Gaussian. The addition of these three Gaussian 
profiles to the broken power law model improves 
the fit, but significant deviations remain in the 
2-6 keV range. These residuals can be improved 
with an extremely broad Gaussian with a peak at 
E = 4.9 keV, a width <r=1.2 keV, and an equiva- 
lent width EW=2.4 keV. (FigureEJ). Even though 
this yields an acceptable fit, the line widths and 
equivalent wi dth seem to be unusua lly extreme for 
a broad line ( Guainazzi et alJl200(J l. although the 
presence of such a line was suggested in earlier ob- 



servations (Gondoin et al. 2002: Longinotti et al 



2007al ). This additional component is not neces- 
sary when a partial covering model is applied to 
the data. 

The narrow 6.4 keV Fe line is present in the 
high- and low-flux spectra. The fluxes in the Fe 
Ka 6.4 keV line in the 2006 and 2007 XMM- 
Newton observations are (1.83±0.33) x 10~ 16 W 
voT 2 and (1.15±0.45)xl0 -16 W m^ 2 , respectively. 
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This result suggests that the flux in the Fe Ka 6.4 
line has become fainter in the low X-ray flux state 
compared to the high state. The equivalent widths 
(EW) however, are different due to the much lower 
continuum flux during the 2007 observation. We 
found EW = 100±20 eV and 200±80 eV dur- 
ing the 2006 and 2007 observations, respectively. 
While the 6.4 keV Fe Ka line is clearly present in 
both spectra, the 7.0 keV Fe XXVI Lya line that 
was found in the 2006 XMM-Newton observation 



(jO'Neill et al.ll2007t ) is absent during the low state 
in 2007. 

Figure[6] displays the difference spectrum be- 
tween the 2006 high-state and 2007 low-state fitted 
in the 2-5 keV energy range with an Galactic ab- 
sorbed power law (ox=l-3) and then extrapolated 
to lower and higher energies. This difference spec- 
trum clearly shows that the changes between the 
high and the low state happened primarily at soft 
energies below 1.5 keV. In some ways this plot is 
similar to the high and low-state sp ectral energy 



distribu tions shown in Figure 2 in iGrupe et al 
(|2007al ). Obviously the high and low state spec- 
tra cannot be described simply by a renormaliza- 
tion of a single power law model. Despite the fact 
that most of the flux change between 2006 and 
2007 occurs at lower energies, the soft X-ray spec- 
tral slopes in the 2006 and 2007 XMM-Newton ob- 
servations are similar. From a broken power law 
model fit the soft X-ray spectral slopes of the 2006 
high and 2007 low states are ax=l-73±0.01 and 
ax=l-89±0.04, respectively. At higher energies, 
however, the energy spectral indices are signifi- 
cantly different. During the 2006 XMM-Newton 
observation the slope was about ax.=1.0, which is 



quite typical for a NLS1 (e.g. iLeighlv |[l999ah . 
The spectral slope in the 2007 XMM-Newton ob- 
servation is very hard with ax=— 0.1. We investi- 
gate two models that can explain such a spectral 
shape: a partial covering absorber and reflection 
from an ionized accretion disk. 

3.2.1. The Partial Covering Absorber Scenario 

Partially covering a single power law contin- 
uum significantly improves the fit to the low-flux 
2007 spectrum compared to the broken power law 
model. The models still requires three Gaussian 
profiles to fit the low energy emission features 
and the narrow 6.4 keV emission line (x 2 /V = 
380/299). The broad Gaussian profile around 5 



keV used in the broken power law model is not re- 
quired. In this scenario the intrinsic power law has 
an energy spectral slope ax=l-78. The absorber 
has a column density of JVh = 15-1 x 10 22 cm~ 2 
and a covering fraction / c =0.94. These partial 
covering parameters are comparable to the ones 
found by S wift during the 2007 May and June ob- 



servations (jGrupe et al. I l2007ah . The partial cov- 



ering parameters for the 2006 observation are then 
AT H pc = 5.5 ±0.2 x 10 22 cm" 2 and / c =0.32±0.01. 

Even though a comparison of the 2007 XMM- 
Newton data with the 2007 May and June Swift 
observations suggests some variability in the col- 
umn density, the spectral parameters are consis- 
tent during these low-flux periods. Notably, all 
spectra fitted by the partial covering absorber 
model only require a single underlying power law 
spectrum. The spectral slope during all these low- 
state observations remains fairly constant within 
uncertainties at ax=l-8. The column densities 
were always within a factor of two. The cover- 
ing fraction at all four epochs was / c =0.94. The 
column density was highest during the 2007 May 



17 Swift observation with iVn.pc = 2 x 10 cm 
and lowest on May 2 5 with iVn,pc = 1 x 10 23 cm 



(jGrupe et al. 2007b). The column densities dur- 



ing the 2007 June Swift and 2007 July 10 XMM- 
Newton observations are in between those values 
with -/Vh, P c = 1.5 x 10 23 cm -2 . The covering frac- 
tion during all these four observation was / c =0.94. 

One way to test whether the partial covering 
absorber model is self-consistent, we fitted the 
2006 and 2007 pn data simultaneously in XSPEC. 
This test can show if the same underlying con- 
tinuum spectrum can be applied to the data only 
affected by the absorber column density and cov- 
ering fraction. The first test was to use a model 
with two power laws (replacing the broken power 
law) and one partial covering absorber. While the 
power law parameters were tied, the partial cov- 
ering absorber parameters were left free to vary, 
first only the covering fraction and second also 
the absorber column density. In none of these 
cases could we obtain an acceptable fit to the spec- 
tra. This result suggests that high and low states 
have different underlying continua. The situation 
changes when adding an additional partial cover- 
ing absorber to the model (XSPEC model wabs * 
zpefabs * zpefabs * (powl + powlj). This model sig- 
nificantly improved the fit from a reduced x 2 — 2-67 
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using the one partial covering absorber to x 2 =l-58 
with the two partial covering absorbers. One of 
the partial covering absorber components, how- 
ever, requires a column density of 3xl0 24 cm~ 2 . 
The covering fractions, which were left as a free 
parameter, of this component are / c =0.2 and 0.8 
for the high and low states, respectively. 

As mentioned by IGrupe et al.1 ( 2007bl ). when 
correcting the fluxes during the low-state obser- 
vations for Galactic and intrinsic absorption, the 
intrinsic fluxes are similar and the intrinsic long- 
term X-ray flux amplitude var ies by a factor of 
about 4-6 (|Grupe et alJ l2007bj)_ which is similar 



to other NLS1 galaxies (e.g. IGrupe et al.l 12001) 



The unabsorbed fluxes in the 0.2-2.0 keV and 2.0- 
10.0 keV band during the 2006 XMM-Newton ob- 
servation were 6.7 and 6.4 xl0~ 14 W m~ 2 , re- 
spectively, which is very similar to the low-state 
XMM-Newton observation on 2007 July 10, which 
were 6.9 and 5.4 xl0~ 14 W m~ 2 , respectively. 
Note that during the extreme low-state on 2007 
May 17 observed by Swift, the fluxes were 4.1 and 
3.2 xl0~ 14 W m~ 2 , respectively. These values 
do not suggest significant intrinsic flux variability. 
Note that the absorption-corrected rest-frame 0.2- 
2.0 keV flux during the ROSAT All-Sky Survey 
fRASS: IVoees et al. Ill999h was 4xl0~ 14 W m" 2 
(Table[3]) , very similar to the flux found during the 
2007 May 17 Swift observations. 

3.2.2. A reflection interpretation for Mkn 335 

The 2007 low-flux, broadband X-ray spec- 
trum of Mrk 335 can be well fitted with a single 
power law continuum and blurred ionized reflec- 
tion (x 2 =1.1). It has been shown that in some 
NLSls the same low-flux model can be applied to 
the high-flux state by simply re-normalizing the 
power law component (e.g. Fabian et al. 2004; 
Gallo et al. 2007a). However, this is not the 
case for Mrk 335 as re-normalizing the power law 
results in a poor fit to the 2006 high-flux state. 
Allowing both the power law and reflection com- 
ponents to vary freely cannot simultaneously de- 
scribe the high-energy (E = 2 — 10 keV) and low- 
energy (E < 2 keV) spectra of Mrk 335 during the 
high-flux state. Similar difficulties were encoun- 
tered by O'Neill et al. (2007) with the same data 
set and demonstrate the need for an additional 
continuum component during the high state. 

One possibility proposed to describe the flaring 



in I Zw 1 (jGallo et al.ll2007allbh was that different 
parts of the accretion disk are illuminated by dif- 
ferent continuum spectra. To model this scenario 
we introduce a second power law with (Model A) 
and without (Model B) associated reflection. In 
Model A we assume that a steep power law (with 
photon index Ti) illuminates the inner accretion 
disc from r m to r^r, while a flatter power law (]?2) 
is associated with outer disc reflection from r^ to 
r out . In Model B, no reflection is associated with 
the steep power law, which may be explained if 
that component is beamed away from the disc; on 
the other hand, in Model B we assume a broken 
power law emissivity profile for the single reflector 
with an inner emissivity index q ln , breaking to q QU t 
at rbr- The two different continuum slopes may in 
fact be due to different mechanisms (e.g. bulk mo- 
tion Comptonization and standard thermal Comp- 
tonization, see e.g. Niedzwiecki & Zdziarski 2006) 
In addition, a MEKAL hot diffuse gas compo- 
nent with Fe abundance tied to that of the reflec- 
tor^) is included to account for any distant emis- 
sion due to thermal plasma which should account 
for the ionized Fe L lines tentatively detected in 
the RGS during the lower-flux state. The param- 
eters of this component are linked in all observa- 
tions. As with any other model discussed above, 
we also include a narrow Fe Ka emission line at 
~ 6.4 keV and a ~ 7 keV line, which is however sig- 
nificantly detected only in the high flux state 2006 
XMM-Newton observation. A narrow Ne line with 
constant intensity at ~ 0.9 keV is also included in 
all fits. 

Both Model A and B describe the various flux 
states well (Figure[7]). See Tab. [6] for fit parame- 
ters. In both models the difference between the 
high- and low-flux state is primarily attributed 
to the significantly diminished flux of the steeper 
power law component. While the two models dif- 
fer in their interpretation of the second power law 
(e.g. beamed versus un-beamed continuum emis- 
sion) , these are issues that can be tested with fur- 
ther analysis (Gallo et al. in prep.). Both models 
are successful in describing the spectra, but they 
both require additional components to explain the 
soft excess. 

3.2.3. Hard X-ray BAT spectrum 

Mrk 335 is detected in the 22 month Swift Burst 
Alert Telescope (BAT; ?) AGN survey (?). The 
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BAT spectrum can be fitted by a single power 
law with ax=l-77io]e4- This is a similar slope 
as found in the Swift XRT and XMM-Newton 
low-state data fitted with a single power law and 
partial covering absorber. Figure[5] displays the 
XMM-Newton pn and Swift BAT spectra fitted si- 
multaneously in XSPEC with a power law with 
a partial covering absorber with ax=l-85±0.03, 
AT HiPC = (16.4 ± 1.0) x 10 22 cm" 2 , and a covering 
fraction of / c =0.94±0.01 (x 2 /z/=582/315). The 
normalizations of the power law models were left 
free. Note that the BAT spectrum is the average 
spectrum over 22 months, so it is not simultane- 
ous with the XMM-Newton low state spectrum. 
The observed 20-40 keV and 20-100 keV fluxes 
are 7.6xl0~ 15 W uT 2 and 1.33xl0~ 14 W m" 2 , 
which is in agreement with the meas urements from 



the 2006 June Suzaku observation (Larsson et al 



20081 ). We do not know from the BAT data if the 
flux in the high-energy bands was lower during the 
low-state in 2007 May to July. However, the AGN 
in the 22 month survey do not show strong vari- 
ability compared with what is typically observed 
at lower energies (J. Tueller, priv comm, 2008). 

3.2.4- Soft X-ray emission lines in the RGS 

The X-ray slopes obtained from single power 
law fits to the EPIC pn (0.2-2.0 keV) and RGS 
(0.35-2.2 keV) spectra are consistent (ax=l-7) for 
the 2006 and 2007 observations. 

The low-flux, broad-band X-ray spectra of Mrk 
335 in Figured] shows significant residuals below 
2 keV when fitted with a single power law. The 
moderate-resolution spectrum suggests the pres- 
ence of an ionized gas either in absorption or emis- 
sion, or both. Figure[9] shows the higher resolu- 
tion RGS spectrum of Mrk 335 fitted with a sin- 
gle power law model. The residuals from this fit 
clearly show the presence of strong emission lines 
around 13, 19, 22, and 25 A. 

The combined RGS 1 and RGS 2 spectrum is 
shown in the lower four panels of Figure[9] in the 
wavelengths bands: 7-14A, 13-2LA, 20-28A, and 
27-34A. Possible identification of the most likely 
features are reported. Most of these lines seem to 
be highly ionized Fe lines including some Fe group 
elements like Mn or Ni. We could also identify 
Ne IX at 13 A. This line, however, is not as strong 
and prominent as was found in the spectrum of the 
partial covering candidate and highly optically po- 



larized NLS1 Mrkl239 (|Grupe et al.ll2004bf) . The 
most prominent group of lines that we have iden- 
tified is the Helium Ka complex of O VII emission 
lines around 22 A. These lines can potentially be 
used for line diagnostics to measure properties of 
the line emitting plasma (§4.2). 

3.3. UV/Optical Data 

3.3.1. UVOT and OM photometry 

The UV/optical photometric measurements of 
Mrk 335 made with the XMM-Newton OM and 
Swift UVOT are listed in TableH The values 
given in the table are c orrected for Galactic red- 
dening (£; B -v=0.035; ISchlegel et all Il998f ) and 



the OM magnitudes were adjusted, as described in 
§2, to be comparable with the Swift UVOT magni- 
tudes. The data suggest that Mrk 335 was slightly 
brighter in the Optical/UV in 2000 than in 2007. 
However, the June 2007 Swift and XMM-Newton 
measurements suggest some rebrightening. The 
origin of the relatively large deviation in the UV 
Wl magnitude during the 2007 XMM-Newton ob- 
servation is not obvious. The 5 reference stars 
listed in Table[2] show a large scatter between the 
XMM-Newton and Swift observations ranging be- 
tween 0.47 to 0.96 mag. However, they all suggest 
that the UV Wl magnitudes between the OM and 
UVOT have to be adjusted significantly. 

As mentioned in iGrupe et al. ( 2007a ). the op- 
tical to X-ray spectral slope a a E of Mrk355 
was between 1.9 and 1.65 during the X-ray flux 
low-state. Ac c ording o the relation given in 



Strateva et al. I (|2005l ) a a ox =1.3 is expected, 
which is e xactly what has been seen durin g the 
high-state (jGrupe et al.ll2007at iGallo 1120061) . 



3.3.2. Optical Spectroscopy 

The upper panel of FigureHU] displays the opti- 
cal spectrum of Mrk 335 taken on 2007 September 
07 at the Xinglong Observatory. The spectrum ap- 
pears nearly identical compared to the spectrum 
taken 8 years earlier at the 1.52m telescope at 
ESO La Silla (jGrupe et al.ll2004al) . The contin- 



uum slope and emission line fluxes are consistent 
within the errors. In Table[7] line measurements 
are given for features between the He II4686 and 



2 The X-ray loudness is defined by Tananbaum et al. I il979f) 



as q ox =-0.384 log(/ 2kcV // 2500 ^). 
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Fe X6375 emission lines. The line measurements 
in the H/3 and [OIII]5007 l ines are, within the er- 
rors, the same as given bv lGrupe et al. ( 2004af ). 

Mrk 335 displays strong coronal Fe lines of [Fe 
VIIJ5159, [Fe XIV]5303, [Fe VII]5721, 6087 and 
[Fe X]6375 (Figure fTOj) . Even though coronal Fe 
line are found of t en in Seyfert 1 galaxies (e.g. 
lErkens et a3ll997l : iMullanev fe Wardl 120071) . see- 
ing these optical lines so strong, as in Mrk 335, 
is rather rare. Similarly strong coronal Fe line 
has been found in the NLS1 galaxies Mrkl239 and 
Ark 564 dErkens et allll997h. and in the Seyfert 1.5 



MCG-06-15-30 ([Reynolds et all 119971) . Unusual 



for a NLS1, Mrk 335 also exhibits strong HeII4686 
and HeI5876 lines in the present spectrum, which 
was also detected previously (e.g., Kassebaum et 
al. 1997). T ypically, this feature is faint in NLS1 
galaxies (e.g. lGrupe et al.ll2004al ). 



4. Discussion 

We reported the XMM-Newton observations of 
the NLS1 Mrk 335 obtained in 2007 July 10, con- 
firming the low-state of Mrk 335 fou nd from Swift 



observ ations in 2007 May/ June (jGrupe et al 
l2007bh . Mrk 335 is remarkable in the sense of 
being still bright enough in its low-state to allow 
a detailed spectral analysis. 

The X-ray spectra of a substantial fraction of 
(NL)S1 galaxies are complex (e.g., Gallo 2006), 
and different types of models have been considered 
to explain them, including reflection (e.g., Ballan- 
tyne et al. 2001, Crummy et al. 2006), partial 
covering (e.g., Tanaka et al. 2005) and relativis- 
tically smeared ionized absorption (e.g., Done et 
al. 2007). Variability induced by (cold) absorp- 
tion is quite common in AGN, and most frequently 
observed in Seyfert 2 galaxies and intermediate 
type Seyferts. An extreme example is the galaxy 
NGC 1365 which shows highly variable cold ab- 
sorption along the line of sight on the time scale 
of days (Risaliti et al. 2007). 

Partial covering models have been frequently 
applied to the X-ray spectra of AGN and NLS1 
galaxies in particular (e.g., Grupe et al. 2004b, c, 
Gallo et al. 2004, Pounds et al. 2004, Tanaka 
et al. 2005), and also to Galactic X-ray bina- 
ries (e.g., Brandt et al. 1996, Ding et al. 2006). 
However, origin and confinement of the required 
dense blobs of gas in AGN partial coverers are 



not well understood (e.g., Kuncic et al. 1997). 
Reflection models generally provide similarly suc- 
cessful X-ray spectral fits. Apart from spectral 
modeling, temporal variability is an important 
discriminant between these two types of models. 
The d eep low state of M rk 335 discovered with 
Swift (|Grupe et al.ll2007rJ) and still present dur- 
ing our follow-up XMM-Newton TOO observation, 
enables us to do a detailed comparison of partial 
covering and reflection models. 

4.1. X-ray continuum 

Before we compare the models themselves we 
compare Mrk 335 with several potentially similar 
NLS1 galaxies. The low-state X-ray spectrum of 
Mrk 335 looks in many ways like t he one seen dur 



ing t he low- state in NGC 4051 (jGuainazzi et al 



19981 ). While iPounds et alj (120041) argue for a par - 
tial covering absorber scenario, Ponti ct al. (2006) 
argue for a reflection dominated spectrum. While 
a fit to the mean low-state spectrum of NGC 4051 
using a partial covering or a reflectio n model re 



suits s tatistically in very similar fits, iPonti et al 



(2006) could show that a reflection model can re- 
produce their findings in the rms spectra and flux- 
flux plots better than a partial covering absorber 
model. Due to the relatively short length of our 
XMM-Newton ToO observation the low-state data 
of Mrk 335 are not of sufficient quality to do such 
an analysis. 

Mrk 335 also displays many similarities to 
the NLS1 galaxy Mrkl239: a) an usual X-ray 
contin uum spectrum, b) strong X-ray emission 
lines ( Grupe et al. 2004rJ). and c) stron g opti- 
cal Fe coronal lines (jErkens et al 1997 ). The 
X-ray spectrum of Mrkl239 was interpreted by 
Grupe et al. (|2004bJ) with a partial covering ab- 
sorber. While Mrkl239 is highly pola rized with 



more than 4% at optical wavelengths ([Goodrich 



19891: iGrupe et al.ll2004bT) . Mrk 335 do es not show 



any significant degree of polarization ( Berriman 



19891 : iBerriman et al.|[l99a ISmith et al 



2002) 



The column density of the partial coverer we 
infer from our spectral fits is similar to those 
often observed in broad absorpti on line (BAL) 



quasars (e.g. IGrupe et al.1 120031 ). In this con- 



text, it is interesting to point out the potential 
similarity of Mrk 335 with the highly variable 
NLS1 galaxy WPVS007 which dropped by a fac- 
tor of more than 400 between the ROSAT All Sky- 
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Surve y observation and follow-u p X-ray observa- 
tions (jGruoe et al.lll995l l2007ah . Neither NLS1 
galaxy shows any detectable changes in their opti- 
cal spectra even though their X-ray fluxes changed 
substantially. While Mrk 335 lacks UV spec- 
troscopy during the low state, the UV spectrum of 
WPVS 007 showed the emergence of strong BAL 
absorption features (jLeighly et al. in prep. 
Recent Swift observations of WPVS 007 suggest 
that the X-ray spectrum is described by a partial 
covering absorber (Grupe et al. 2008, in prep). 

An interesting connection can also be drawn to 
the X -ray weak quasar PG 2112+059 (jSchartel et al 
l2007h . Also here the X-ray spectra can be fitted 
by partial covering absorber and reflection mod- 
els. Although the fits are statistically similar (like 
what we saw in Mkn 335) the authors prefer the 
blurred X-ray ionized reflection model, because 
the absorber model yields unreasonable high col- 
umn densities and ionization parameters. In gen- 
eral, studying bright AGN in low-state like Mrk 
335 gives us access to well-exposed X-ray spectra 
of X-ray weak AGN which appear to be signif- 
icantly fainter in X-rays th an what is expecte d 



from optical luminosities (jBrandt et al. 2000). 



Typically these AGN are too faint in X-rays to 
obtain well-exposed X-ray spectra. Note that the 
situation in Mrk 335 is different compared with 
the X-ray weak qu asar PHL 1811 whic h is intrin- 
sically X-ray weak (jLeighly et al.ll2007l ). 



Spectral Fits. The spectral fits to the X-ray 
spectra of Mrk 335 involving a partial coverer are 
relatively simple, in the sense of only involving 
two free parameters to describe the partial cov- 
erer itself. All low-state spectra of Mrk 335 can 
be fitted by the same underlying single power law 
with Q!x=1.8. However, the XMM-Newton high 
state data cannot be fitted by the same under- 
lying power law. The 2006 XMM-Newton data 
do require an underl ying broken power law model 
(jGrupe et al ] l2007bh . We performed several tests 
with double power law with one and two partial 
covering absorbers. We found that we can fit the 
high and the low states with the same underly- 
ing continuum model, but the low state spectrum 
then requires a very high column density of several 
times 10 24 cm -2 to describe the spectrum at ener- 
gies > 2 keV. On the other hand, a constant flux 
of the underlying continuum may explain why the 



optical spectrum has not changed at all between 
1999 and 2007. 

The steep X-ray spectral slope ax=l-8 contin- 
ues even into the Swift BAT hard X-ray energy 
range. While X-ray spectral slopes of ax=l-8 are 
quite common for NLSls in the 0.2-2.0 keV range 
(|Boller et al. Ill996l : iGrupe et al.ll200lh . this steep 
X-ray spectral slope is unusual at hard X-ray en- 
ergies. In the BAT survey (?) no AGN has been 
found with an X-ray spectral slope as steep as Mrk 
335. Note that also in the 2-10 keV range an X-ray 
spectr al slope ax=l-8 is quite unusual ( ^Leighly 



1999ri ). As suggested by e.g. ? steep hard X-ray 
spectral slopes can be caused by thermal Comp- 
tonization of UV photons of the hot accretion disk 



X-ray and UV variability. The short-time 
variability found during the low-state observa- 
tion shows that the bulk of the soft X-ray emis- 
sion is not fro m an extend e d scatt ering region as 

" (12004 ' 



for Mrkl239 
for Mrk6. These re- 



su ggested by IGrupe et al 



or 



Immler et al. 



(I2OO3I) 



gions are typically too large to show variability 
on time scales of hours. Note that the mass 
of the central black hole i n Mrk 335 is of the 
order of 1.5 - 3 x 10 7 M^ dPeterson et all [2004: 



Vestergaard fc Peterson I |2006J) , which is on the 
high-mass end of the bl ack hole mass distribu- 
tion of NLS1 gala xies (|Grupe fc Mathurl 12004 
Watson et alll2007l ). 

The variability detected during the low-state 
observation implies a source size of a few light 
hours (see Figure[2|) which in turn requires that 
the partial coverer has to consist of very compact, 
dense clumps of gas. 

Perhaps in favor of the partial covering inter- 
pretation is the variability in the UV/optical band 
(even though small). In the two bluest filters, W2 
and M2, Mrk 335 varies stronger by 0.2 mag than 
in the optical V and B filters. Since the UV mag- 
nitudes seem to follow the X-ray count rates, this 
might indicate that the variability in the UV and 
in X-rays is caused by the same mechanism. How- 
ever, longer-term monitoring of the UV and X-ray 
variability is required to put this statement on a 
more firm basis. 
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Models of partial coverers. In summary, the 
most successful phenomena-logical partial coverer 
model consists of very compact blobs of gas cov- 
ering the variable X-ray source and is mostly lo- 
cated along our line of sight in order to ensure non- 
variable BLR line emission and in order to be con- 
sistent with the absence of a detectable strongly 
variable 6.4 keV iron line from the partial cov- 
erer itself. The requirement of the compactness of 
the clouds, in combination with their high column 
density, emphasizes well-known challenges to the 
partial coverer scenario including questions like: 
what is the origin of the blobs, and how are they 
confined (e.g., Guilbert & Rees 1988, Celotti et 
al. 1992, Kuncic et al. 1997) ? In the model 
of Abrassart & Czerny (2000), fragmentation of 
the inner accretion disk leads to the formation of 
compact clouds which can partially obscure the 
central X-ray source, confined perhaps by strong 
magnetic fields (Kuncic et al. 1997). These clouds 
would reprocess some radiation into the optical 
and UV band (Celotti et al. 1992, Kunic et al. 
1997, Abrassart & Czerny 2000). If their emission 
does contribute to the total UV emission, then the 
variability in the X-ray band would be higher than 
that in the UV band, as observed. 

Reflection models. Blurred reflection models 
invoking a single continuum source and reflecting 
medium have been sh own to fit the 2000 high- 
flux state observation ( Crummy et alJ lioOrjf ) and 
the 2007 low-flux state observation (Section 3.2.2). 
However, there are a number of discrepancies be- 
tween parameter values (e.g. inclination and iron 
abundance) that indicate these models arc not 
consistent. Likewise, the high signal-to- noise 2007 
spectrum of Mrk 335 cannot be modeled self- 
consistently in the high- and low-energy ranges 
(O'Neill et al. 2007; Section 3.2.2). 

The fact that the low-flux spectra of Mkn 335 
cannot be fit with the same reflection model 
used for the higher flux states suggests that the 
most straightforward applic ation of light bending 



does not necessarily app ly (jMiniutti et alJ 1200 



Miniutti fc Fabian lliool . That is, the low-flux 
state is not solely due to a single continuum source 
whose intensity appears diminished because of 
its proximity to the black hole while the reflec- 
tion is kept constant. One possibility is that the 
entire disc does not "see" the same continuum 



source. Such a scenario was proposed by Gallo et 
al. (2007a, b) to describe the behavior of another 
NLS1, I Zw 1. Here we considered this possibility 
for Mkn 335. 

Two models are suggested, both rendering 
equally good fits. Both scenarios incorporate two 
power law components, each one illuminating dif- 
ferent parts of the accretion disc. One component 
could be associated with a jet producing either 
beamed (Model B) emission with no reflection 
component, or un-beamed (Model A, see Ghis- 
cllini et al. 2004) emission that preferentially 
illuminates the inner-most region of the disc. A 
jet seen in Mrk 335 is plausible. Mrk 335 has a 
clear detection of a compact source at 4.8 GHz 
2006f) . The second power law is as- 



([Leipski et al 



sociated with the traditional idea of the corona, 
which blankets the disc to larger radii. These, of 
course, are not exclusive interpretations. The two 
different continuum slopes may in fact be due to 
different mechanisms (e.g. bulk motion Comp- 
tonization and standard thermal Comptonization, 
see e.g. Nicdzwiecki & Zdziarski 2006). 

Statistically, both models fit the data very well 
(xt — 1-06 — 1-10), but more importantly in a 
self-consistent manner. The primary difference be- 
tween the high- and low-flux states is the dimin- 
ished flux received by the observer from the inner 
power law component (i.e. the one perhaps asso- 
ciated with a jet). 

According to our models A and B, there are sit- 
uations (especially in the low state) when the X- 
ray spectrum of Mrk 335 is reflection-dominated. 
If so, the observed flux cannot be used in a sim- 
ple way to infer the intrinsic X-ray luminosity of 
the AGN. In general, the power law component 
contributes to the observed flux at infinity, and it 
also illuminates the accretion disc giving rise to 
the associated reflection component. The intrinsic 
luminosity of that power law must take into ac- 
count not only the flux observed at infinity, but 
also that needed to produce the reflection compo- 
nent. If i D bs is the power law luminosity observed 
at infinity and R is the reflection fraction of the as- 
sociated reflection component, the intrinsic power 
law luminosity is L mt — L \> S (R + l)/2, which, in 
the isotropic case (R = 1), simply reduces to I/ bs- 
Applying the above arguments to model A, we es- 
timate that the intrinsic luminosity of all power 
law components is ~ 2.3 x 10 43 erg s _1 cm" 2 in the 
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high-flux state and >= 1.9 x 10 43 erg s _1 cm' 2 in 
the low-flux state (all luminosities are unabsorbed 
and given in the 2-10 keV band). In other words, 
the intrinsic luminosity of the AGN would be simi- 
lar in the two states despite the dramatic observed 
flux variation. 

In the case of model B, we have one power 
law with its reflection (for which the same argu- 
ments as above can be applied) and one power law 
with no associated reflection. The latter could 
be due to a jetted component beaming the ra- 
diation away from the disc (so that there is no 
reflection). Hence, we can estimate its intrin- 
sic luminosity by assuming that the radiation is 
beamed in the half-sphere away from the disc 
(Lint <= Lobs/2). For model B, considering the 
contributions of both the jetted and the non-jetted 
components, we then obtain intrinsic luminosities 
in the range (6 — 6.8) x 10 43 erg s _1 cm -2 in the 
high-flux state and ~ 2.4 x 10 43 erg s _1 cm~ 2 in 
the low-flux state (the jetted contribution is neg- 
ligible in the low- flux state). 



Partial covering absorber vs. reflection 
models. Statistically the reflection models yield 
better fits compared with the partial covering ab- 
sorber model (reduced x 2 =1.10 vs. 1.27, respec- 
tively). Both, the partial covering absorber model 
and the reflection models require underlying bro- 
ken power law or double power law model spectra 
to fit the high-states of Mrk 335. For both models 
the same underlying continua can be used in the 
high and low states just modified by either absorp- 
tion or reflection. The parameters however, used 
for the partial covering absorber model then re- 
quire that the absorption column density of one of 
the partial covering absorber components is large 
on the order of several times 10 24 cm~ 2 . Even 
though this is not impossible, it requires extremely 
dense compact clouds to pass through the line of 
sight. 

Reflection models have been used in the past 
to describe t he X-ray spectra of may AGN (e.g. 
iFabianl 12007ft and have also been applied suc- 
cessfully to describe the high-state X-ray spec- 



possibility - that the spectrum is dominated in 
high-state by reflection, and this reflection spec- 
trum is modified by partial covering absorption. 
A scenario like t his h as been suggested e.g. by 
Chevallier et al. I (|2006h and ?. In some ways the 
situation in Mrk 335 is similar to that of PG 
2112+059 where the absorber model also yields 
very high column densities which seems to fa- 
vor the reflection model. On the other hand the 
partial covering absorber model is rather simple 
and does not require many parameters to fit the 
data reasonably well. At the end one has to con- 
clude that both models are valid to explain the 
spectral shape seen in the low state spectrum of 
Mrk 335. Nevertheless with the low-state XMM- 
Newton observation presented here we were able 
to put strong constrains on both the partial cov- 
ering and reflection model parameters. In order 
to really distinguish between both models either a 
very long XMM-Newton observation in low state 
is needed to search for the deviation between both 
models at intermediate energies between 2-6 keV 
or it is necessary to obtain a spectrum at energies 
> 10 keV using e.g. Suzaku. 

4.2. X-ray emission lines 

Several low-energy (E < 2 keV) X-ray emis- 
sion lines are detected in the RGS spectra of 
Mrk 335. Soft X-ray emission lines are typi- 
cally found in obscured Seyf er t 2 type AGN (e.g. 



Guainazzi & Bianchi I 2007: Kinkhabwala et al 



2002) and are only seen in Seyfert Is when the 
AGN enters a low-flux state. Consequently, low- 
energy emission lines have only been detected in 
a few Seyfe rt Is like: Mr k590, NGC 4051, and 



NGC 5548 (ILonginotti et alj 2007ct Pounds et al . 
2004t ISteenbrugge et al.ll2003L respectively). This 



tra of Mkn 335 (e.g. iBallantvne et al 



Crummy et al.l2006tlLonginotti et al.l2007at 



2001 



20071) and we have shown here that they also con- 
sistently produce the low-state. Further, we men- 
tion in passing that we cannot exclude a third 



O'Neill et al 



is the first time such lines have been observed in 
Mrk 335. Even though these lines are expected in 
Seyfert Is when the continuum is suppressed, as 
in the low state 2007 XMM-Newton observation, 
they are very prominent even with just a 22 ks 
RGS exposure. 

Most notable in the RGS spectrum of Mrk 335 
is the O VII emission line complex around 22A. 
The forbidden (f) line at 22.1012 A and the inter- 
combination (i) line at 21.8044 A are secure iden- 
tifications. These lines have also been found in the 



low-s t ate RGS spectra o f NGC 4051 ([Pounds et al 
2004t IPonti et al.l l2006h . The resonance (r) line 
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appears to be weak and slightly red-shifted with 
respect to its rest frame wavelength of 21.602 A, 
which could be due to resonance scattering and 
absorption. Although the line intensities are con- 
siderably uncertain, the presence of the O VII line 
complex suggests coronal emission from collision- 
ally ionized plasma. In principle the lines from 
this complex can yield reasonably precise infor- 
mation on the temperature, densities, and ioniza- 
tion state of the source. However, owing to the 
uncertainties we conclude that it is not useful at 
this time to attempt a detailed line ratios anal- 
ysis such as t he G = (i+f)/r ratio or the R = 



f/i ratio (e.g. iPradhan fe Shull Ill98lh . Nor is it 



possible to deduce the existence of additional fea- 
tures widely se en in warm absorber spectra (e.g., 
Lee et al.ll200ll ). such as the OVI KLL absorption 
features at (22.05,21.87) A lying between the f 
and i lines (|Pradhan I fcoOOl : IPradhan etaD 120031 ) . 
In addition to O VII, we have identified, to vary- 
ing degree of certainty, a number of iron lines 
from coronal ions such as Fe XVII to Fe XXV. 
The iron lines complement the independent ob- 
servations of the UV/Optical lines from coronal 
iron ions in lower ionization stages discussed in 
the next section. Note, that the RGS data have 
ndepe ndently been analyzed by lLonginotti et "al 



( 20081 ) who present line ratios. They did not find 
a unique solution of the plasma properties either, 
but concluded that the location of the X-ray line 
emitting gas is within 0.01 pc. 

4.3. Optical emission lines 

The optical emission-line spectrum of Mrk 335 
is remarkably constant between the observations 
taken 8 years ago, and the one taken nearly si- 
multaneously with the X-ray low-state observa- 
tions. Since the broad component of the Balmcr 
lines and the broad component of the Hell line are 
very sensitive to the EUV-to-soft X-ray part of the 
continuum, the observed lack of variability implies 
that the bulk of the BLR still sees the high-state 
continuum. In the context of the partial cover- 
ing model, this can be understood if the absorber 
has a small global covering fraction and is mostly 
located along our line-of-sight. 

The highly ionized coronal Fe lines found in the 
optical spectrum of Mrk 335 are rather strong. 
While coronal Fe lines are present in the opti- 
cal spectra of Seyfert Is, they typically appear 



faint (e.g. vanden Berk et all 2001 ) with typical 
flux ratios to [OIII]A5007A about 0.0 3 to 0.06. In 
compa rison with the AGN sample of lErkens et all 
(|l997t ). only 2 out of 15 AGN (Mrk486 and Akn 
564), have coronal Fe lines stronger than Mrk 335. 
The central wavelengths of the iron lines in Mrk 
335 suggests they are blueshifted with respect to 
the [OIII] lines, as is commonly seen in coronal 
lines (e.g., Erkens et al. 1997, Rodriguez- Ardila 
et al. 2006). The 5A resolution of our spec- 
trum, however, does not allow a precise measure- 
ment of the central wavelength of the lines and 
therefore the shift of the line with respect to the 
rest-frame cannot be determi ned with high accu - 
racy. As generally observed (jErkens et all Il997 ). 
the line widths of the coronal Fe lines are between 
the NLR (FWHM« 500 km s" 1 ) and the BLR 
(FWHM(H/3)=1700 km s^ 1 ). 

5. Conclusions 

We report the XMM-Newton and Swift observa- 
tions of the bright NLS1 galaxy Mrk 335 in a deep, 
and unprecedented, low-flux state. The main find- 
ings are as follow: 

• The XMM-Newton ToO observation of 2007 
July 10 still found Mrk 335 in a very low- 
flux state as previously detected by Swift. 
With this XMM-Newton observation, how- 
ever, the parameters of the X-ray spectral 
models can be well-constrained. Low-state 
observations of Seyfert Is are rare and Mrk 
335 is exceptional in being still bright in its 
low-state. Low states give us access to X-ray 
emission lines that otherwise cannot be seen 
(see below) . The low state of Mrk 335 lasted 
at least two months. 

• While a simple broken power law model re- 
quires a broad red-shifted Fe Gaussian line 
component with an extremely large equiva- 
lent width EW=2.4 keV, this line compo- 
nent is not required either in the partial 
covering absorber nor the reflection mod- 
els. The partial covering model can be fit 
to high-state and low-state data with the 
same underlying double power law model 
continuum when two absorber components 
are used, even though formally a second 
power law component is not required to fit the 
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low-state spectrum. The reflection model 
formally accounts for the two-component na- 
ture of the observed spectrum by multiple 
reflection where different parts of the accre- 
tion disk sec different continua. 

• High-resolution spectroscopy clearly reveals 
the presence of X-ray emission lines in the 
low-energy X-ray spectrum of Mrk 335, 
including transitions of oxygen and iron. 
These lines are only visible in the low- 
state. Preliminary results suggest that 
the emission-line gas is collisionally ionized. 
However, to perform proper line diagnostics 
of the ionized gas to derive gas density and 
temperature requires a longer RGS exposure 
of Mrk 335 in its low-state. 

• The optical spectrum of Mrk 335 taken dur- 
ing the X-ray low-state did not change in 
comparison to a spectrum taken 8 years ear- 
lier. The constancy of the broad Balmer 
lines indicates that the bulk of the broad 
line region was still illuminated by the "high- 
state" continuum. The optical spectrum of 
Mrk 335 displays a variety of strong, highly- 
ionized coronal Fe lines. These might be 
linked to the component emitting the X-ray 
lines detected with RGS. The non- variability 
seen in these lines may indicate that the in- 
trinsic SED has not changed. 

Mrk 335 will continue to be monitored with 
Swift in order to trigger another, much deeper 
X-ray observation if the AGN enters another low 
state. 
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Table 1 

Summary of the 2007 July 10 XMM-Newton observations of Mrk 335 



EPIC pn 


17:09 


21:57 


16779 


EPIC MOS-1 


16:14 


21:56 


20348 


EPIC MOS-2 


16:14 


21:56 


20351 


RGS-1 


15:50 


21:58 


21932 


RGS-2 


15:50 


21:58 


21936 


OM U 


15:36 


16:32 


3000 


OM B 


16:32 


17:10 


2000 


OM UVW1 


17:10 


18:22 


4000 


OM UVM2 


18:22 


19:34 


4000 


OM UVW2 


19:34 


20:46 


4000 


OM UVW2 


20:46 


21:57 


4000 



1 Start and End tii 
2 Obscrving time g 



Table 2 

List of reference stars used to determine the error in the photometry of Mrk 335 and 

THE OFF-SET BETWEEN THE Swift-XJV OT AND XMM-Newton-OM FILTERS. THE MAGNITUDES ARE NOT 

CORRECTED FOR GALACTIC REDDENING. 



XMM 2000 



XMM 2007 



Swift 00 1 1 



Swift 002 1 Swift 003-005 1 



Star 


1 00 


06 


22.8 


+ 20 14 58.0 


V 


14. 


28±0. 


01 








14.29 + 0.01 


14.28 + 0.01 


14.28+0.01 












B 


14. 


94 + 0. 


01 


14. 


90±0 


.01 


14.89+0.01 


14.87+0.01 


14.86 + 0.01 












U 


14. 


88 + 0. 


01 


14. 


88 + 


.01 


14.87+0.01 


14.87+0.01 


14.86 + 0.01 












UV Wl 








15. 


44±0 


.01 


16.27+0.02 


16.32+0.02 


16.33 + 0.02 












UV M2 


17. 


43 + 


16 


17. 


45 + 


.09 


17.90 + 0.06 


17.91 + 0.05 


17.94 + 0.05 












UV W2 








17. 


68 + 


.20 


17.73+0.04 


17.78+0.03 


17.81 + 0.03 


Star 


2 00 


06 


32.0 


+ 20 13 08.3 


V 


15. 


46±0. 


02 








15.48 + 0.02 


15.48 + 0.02 


15.51 + 0.01 












B 


16. 


31 + 


01 


16. 


29±0 


.01 


16.21 + 0.02 


16.20 + 0.01 


16.21 + 0.01 












U 


16. 


53 + 0. 


02 


16. 


54 + 


.02 


16.55 + 0.03 


16.57+0.02 


16.55 + 0.02 












UV Wl 








17. 


26 + 


.03 


18.30+0.07 


18.20+0.05 


18.17 + 0.05 












UV M2 














>20.03 


>20.35 


>20.42 












UV W2 














19.76 + 0.12 


19.76+0.10 


19.91 + 0.10 


Star 


3 00 


06 


17.1 


+ 20 14 08.4 


V 


15. 


40±0. 


02 








15.46 + 0.02 


15.43+0.02 


15.42 + 0.01 












B 


16. 


21±0. 


01 


16. 


19 + 


.01 


16.13 + 0.02 


16.13 + 0.01 


16.09 + 0.01 












U 


16. 


27±0. 


02 


16. 


26 + 


.02 


16.26+0.02 


16.26+0.02 


16.23+0.02 












UV Wl 








16. 


88 + 


.03 


17.69 + 0.05 


17.83 + 0.04 


17.74 + 0.0 1 












UV M2 














19.85 + 0.19 


19.82 + 0.15 


19.69 + 0.13 












UV W2 














19.27+0.09 


19.34+0.07 


19.45 + 0.07 


Star 


4 00 


06 


18.0 


+ 20 13 17.2 


V 


15. 


04+0. 


01 








15.04+0.02 


15.07 + 0.01 


15.07+0.01 












B 


15. 


51 + 0. 


01 


15. 


49 + 


.02 


15.46+0.01 


15.45+0.01 


15.44+0.01 












U 


15. 


35 + 0. 


01 


15. 


36 + 


.01 


15.33+0.01 


15.35+0.01 


15.34 + 0.01 












UV Wl 








15. 


70 + 


.02 


16.35+0.02 


16.34+0.02 


16.37+0.02 












UV M2 


16. 


85±0. 


12 


16. 


77+0 


.06 


16.99 + 0.04 


17.02+0.03 


17.09 + 0.03 












UV W2 








17. 


16+0 


.13 


17.35 + 0.03 


17.37 + 0.02 


17.42 + 0.02 


Star 


5 00 


06 


20.1 


+ 20 10 50.2 


V 


14. 


30±0. 


01 








14.31 + 0.01 


14.33+0.01 


14.30+0.01 












B 


15. 


51±0. 


01 


15. 


47 + 


.01 


15.38+0.01 


15.39+0.01 


15.39 + 0.01 












u 


16 


45±0. 


02 


16. 


45 + 


.02 


16.46 + 0.03 


16.41 + 0.02 


16.48+0.02 












UV Wl 








17. 


54+0 


.04 


18.01 + 0.06 


18.12 + 0.05 


18.12 + 0.04 












UV M2 














>20.03 


>20.34 


>20.39 












UV W2 














19.22+0.09 


19.44 + 0.08 


19.52 + 0.07 


1 Thc 


observir 


•g t 




of the 2007 Swift 


observation; 




listed 




Table 


1 inlUrupe 


et alj BOOTH). 
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Table 3 

Summary of the 0.2-2.0 keV and 2.0-10.0 keV fluxes 



Mission Time 1 0.2-2.0 koV flux 2 2.0-10.0 keV flux 2 References 3 



UHURU 


1971/1972 




2.04±0.03 


(1) 


HEAO 1 


1978-07-07 




1.1±0.2 


(1) 


EINSTEIN 


1981-01 




1.46±0.10 


(1) 


EXOSAT 


1983-11-05 


1.4±0.2 


0.23±0,05 


(2) 


EXOSAT 


1984-12-06 


8.3±0.1 


1.54±0.03 


(2) 


GINGA 


1987-12-03 




0.91±0.16 


(3) 


GINGA 


1988-11-29 




2.00±0.10 


(3) 


ROSAT RASS 


1990-07-14 


4.00±0.10 




(4) 


BBXRT 


1990-12-05 




0.84±0.05 


(1) 


ROSAT PSPC 


1991-06-30 


3.40±0.10 




(4) 


ASCA 


1993-12-09 




0.94±0.10 


(5) 


X MAI- Newton 


2000-12-25 


7.90±0.01 


1.58±0.02 


(6) 


XMM-Newton 


2006-01-03 


6.76±0.01 


1.79±0.01 


(6) 


Suzaku 


2006-06-20 




1.43±0.02 


(7) 


Swift 001 


2007-05-17 


0.23±0.03 


0.28±0.04 


(6) 


XMM-Newton 


2007-07-10 


0.41±0.01 


0.31±0.01 


(8) 


1 Observing tim 










2 Thc 0.2-2.0 ar 


d 2.0-10.0 keV flux 


es are corrected for 


Galactic absorption and arc give 


in units of 10- 14 


W m -2 (= lO" 11 


orgs s-1 em" 2 ). 






° References: (1) Turner et al. (1993). (2) Pounds etal. i 19871), (3) 


Nandra h Pounds] 
), (B)IUruoe et al.| 


J200?b|). (7)|Larss 


i; l hi. liuuiij, y^j ) vjeuij^e l-l tLi . i iuuu ^, j_i e i it, 1 1 1 v ijya. 

on et ill . 1 2008 ) . (8) this paper 





Table 4 

Optical and UV photometry of Mrk 335 using the Swift-WOT and XMM-Newton-OM 1 



Filter 


XMM 2000 


XMM 2007 


Swift 001 


Swift 002 


Swift 003 


Swift 004 


Swift 005 


V 
B 

U 

UV Wl 
UV M2 
UV W2 


14.00±0.01 
14.15±0.01 
13.00±0.01 

12.92±0.16 


14.25±0.01 
13.17±0.01 
13.52±0.01 
12.90±0.01 
13.07±0.20 


14.22±0.01 
14. 49 ±0.01 
13.29±0.01 
13.20±0.01 
13.12±0.01 
13.19±0.01 


14.19±0.01 
14.42±0.01 
13.20±0.01 
13.08±0.02 
12.95±0.01 
13.01±0.01 


14.19±0.01 
14.44±0.01 
13.23±0.01 
13.09±0.02 
12.98±0.01 
13.04±0.01 


14.15±0.01 
14.41±0.01 
13.22±0.01 
13.10±0.01 
12.97±0.01 
13.04±0.01 


14.18±0.01 
14.42±0.01 
13.23±0.01 
13.05±0.01 
12.93±0.01 
12.98±0.01 



1 A11 magnitudes are corrected for Galactic reddening (E B _ v =0 . 035 ; Schlegel et al. 1998) and the XMM-Newton OM 
magnitudes are adjusted as described in the text to be comparable with the Swift UVOT data. 
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Table 5 

Spectral analysis of the 2007 and 2006 XMM-Newton spectra of Mrk 335 in the 0.3-12.0 

keV range 



"H.pcf 



XMM-Newton 2007 EPIC pn 


(a) 


1 


63±0.02 


















3158/312 




(b) 


1 


90±0.02 


1 


73±0 


03 


-0.13±0.03 










712/308 






1 


s -+0.03 
° -0.06 
S7 +0.03 
° -0.06 


1 


6S±0 


04 


-0.10±0.04 






6.37±0.05 


0.31±0.13 


505/299 




(d) 


1 


1 


77+° 

-o 


16 
08 


+ 0.09 
+ °- 05 _0.07 






6.42±0.04 


0.19±0.08 


333/296 






1 


TS + 0.02 
'°-0.03 










15 1 + 1 - 
15,1 -0.9 


0.94±0.04 


6.41±0.04 


0.19±0.08 


380/299 


EPIC pn + MOS 




1 


55±0.02 


















5678/549 




(b) 


1 


89±0.02 


1 


73±0 


02 


-0.19±0.03 










1208/547 






1 


7 o + 0.03 
'°-0.02 










14 6 +1 "° 


93+ " 04 


6.41±0.03 


0.26±0.07 


788/538 


RGS1 + 2 




1 


64±0.14 


















424/371 


XMM-Newton 2006 EPIC pn 




1 


62±0.01 


















16770/697 




(b) 


1 


72±0.01 


1 


S4±0 


02 


1.07±0.01 










1717/695 






1 


72±0.01 


1 


79±0 


03 


1.10 = 0.02 






6.41±0.04 


0.10±0.02 


1043/686 






1 


73±0.01 










5.1±0.4 


0.35±0.02 


6.41±0.03 


0.09±0.02 


834/687 


RGS1 + 2 




1 


64±0.01 


















5288/4494 


XMM-Newton 2006 + 2007 EPIC pn 




1 


62±0.01 


















19930/1008 




(b) 


1 


72±0.01 


1 


84±0 


02 


1.06±0.01 










4131/1006 






1 


72±0.01 


1 


67±0 


03 


1.28±0.02 






6.41 (fix) 


0.08/0.16 


1710/997 






1 


70±0.01 










10. 8 7 


0.56/0.91 7 


6.41 (fix) 


0.10/0.30 7 


2642/997 






1 


69±0.01 










10. 5/14. 4 7 


0.56/0.93 7 


6.41 (fix) 


0.10/0.20 7 


2575/996 



-'■Spectral models used are: (a) absorbed power law, (b) absorbed broken power law, (c) absorbed broken power law with two Gaussians to represent the soft X-ray emission 
lines below 1 keV and one Gaussian for the 6.4 keV Pe line, (d) same as (c) but with an additional broad line at 4.9 keV with a cr=1.5 keV, and (c) power law with a partial 
covering absorber with two soft X-ray Gaussians and one for the 6.4 Fe line. For all models the absorption column density of the z = absorber was fixed to the Galactic value 
(3.96X10 20 cm- 2 rTJTcEeV h LockmaiTirTMTD . 

2 The break energy S brcak is given in units of keV. 

3 Absorption column density of the partial covering absorber iV H f in units of 10 22 cm -2 
4 Covcring fraction f c 

5 The line energy E Unc is given in units of keV. 
6 The line equivalent width is given in units of keV. 

Simultaneous fit with a power law with partial covering absorber and Gaussian lines. While the X-ray spectral slope was tied and the line energies and widths were fixed, 
the column densities and covering fractions of the partial covering absorber and the normalizations of the power law continua and the line fluxes were loft as free parameters. 
The values of the column densities, covering fractions, and equivalent widths of the 6.4 keV line are given for the 2006 and 2007 data. 



Table 6 

Summary of the most important parameters 3, for the best-fitting reflection models 
(Model A and B; see text for details) to the high-flux (HF) and low-flux (LF) state 

observations. 



Parana. 


HF state 


LF state 


ticd b 


Units 






Model A 


( X 2 r = 1.06) 







Oil 


i.86^:S 


1.86±££ 




X 


[xlO~ 3 W m" 












(5.86 ±0.05) x lO" 11 


< 3.5 x 10~ 13 






" 2 ] 


= [ergs s 


1 cm 






0.70 ±0.01 


0.70 ±0.01 




X 












L a 2 


(1.93 ± 0.01) x 10~ n 


(6.9 ± 1.3) x 10" 


1 Q 




[xlO~ 3 W irT 


" 2 ] 


= [ergs s" 


1 

cm" 


~ 2 ] 


Tin 


1.6 ±0.1 


1.6 ±0.1 




X 


r c 
'g 










i 


39 ± 1 


39 ± 1 




X 


degrees 










Qi 


< 3 -°-o.6 






X 












<?2 


4.3 ± 0.3 


4.3 ± 0.3 




X 












7"br 


40 -18 


> 1.8 
















£l 


150±?| 


38± 2 






[xlO~ 3 W irT 


■ 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 


F reil 


(3.3 ±0.9) x 10~ 12 


(1.1 ±0.5) x 10" 


12 




[xl0~ 3 W irT 


- 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 


6 


150±i| 


38t 2 






[xl0~ 3 W irT 


- 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 




(3.9 ±0.1) x 10~ 12 


(6.4 ±0.4) x 10- 


12 




[xl0~ 3 W in- 


- 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 


^4fo 1=2 


1.6 ±0.1 


1.6 ±0.1 




X 


solar 












Model B 


( X 2 r = 1.10) 
















on 


1 64 +(u > 4 

1 - D4 -0.02 


1 64 +(u)4 




X 












F ai 


(4.55 ±0.01) x 10~ n 


(9.7 ±5.1) x 10" 


14 




[xlO~ 3 W irT 


- 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 




0.87 ±0.01 


0.74 ±0.01 


















(5.87 ±0.09) x 10~ 12 


(6.27 ±0.63) x 10 


-13 




[xlO~ 3 W irT 


■ 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 


Tin 


< 1.34 


< 1.34 




X 












i 


50 ± 1 


50 ± 1 




X 


degrees 










Qin 


> 4.5 


5.9 ±0.2 
















Qout 


> 4.5 


4.4 ±0.2 


















2.8 ±0.1 


2.8 ±0.2 






r c 
'g 










c 


20 ±1 


24 ±1 






[xl0~ 3 W m" 


- 2 ] 


= [ergs s" 


1 cm" 


■ 2 ] 




(1.23 ±0.01) x 10~ n 


(5.6 ±0.1) x 10" 


12 




[xl0~ 3 W in- 


- 2 ] 


= [ergs s" 


1 cm" 


- 2 ] 


-4fo 1=2 


0.9 ±0.1 


0.9 ±0.1 




X 


solar 











a All disk radii are given in units of the gravitational radius and the outer disk radius is fixed to 400 
r\ rmg . The temperature of the MEKAL component is 1.4 kcV. Other parameters listed are fluxes F, 
the inclination angle i, the Fe abundance (Ap e ) given relative to solar, the ionization parameter (£), 
and the emissivity index (q) 

b This column marks which model parameters have been tied between the high and low-state spectra. 

c Gravitational radius r g = GM/c 2 . 
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Table 7 

Measurements of the optical emission lines 



Line 




EW 1 


*c,ob, ! FWHM 3 


x/[OIII] 4 


He II4686n 




5.0 


4807 


950 


0.19 


He 114686b 




16.9 


4799 


5100 


0.59 


H/3n 




4.0 


4984.2 


560 


0.13 


H/3b 




90.0 


4984.2 


1700 


3.23 


[OIIIJ4959 




8.4 


5086.7 


500 


0.31 


[OIIIJ5007 




28.1 


5136.1 


500 


1.00 


[FeVII]5159 




3.4 


5293.3 


1250 


0.11 


[FeXIV]5303 




2.8 


5441.3 


1350 


0.09 


[FeVII]5721 




4.0 


5866.7 


1250 


0.14 


Ho I 5876 




15.8 


6025.7 


1420 


0.54 


[FeVII]6087 




3.0 


6239.6 


930 


0.09 


[FoX]6375 




3.9 


6535.8 


900 


0.13 


^Thc equivalent widths 


are given in u 


nits of A. 




2 The observ 


ed 




'avelength of e 


I lino is g 




3 Full-width 




half-max 


imum of the 1 


ine given 


in units of 


km s _1 . 












4 Line fluxes 




native to 


[OIII]A5007 


lino flux 


< F [OIII] = 


1.93 X 10~ 16 


W 


m — 2 
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Fig. 1. — Long-term light curves of Mrk 335. The 
upper panel shows the 2-10 keV light curve start- 
ing with the first X-ray detection of Mkn335 by 
UHURU in 1971, the middle panel displays the 
0.2-2.0 keV flux light curve starting at the EX- 
OSAT observation in November 1983. The lower 
panel displays the UV flux light curve measured 
with IUE and HST at 2256A, and XMM-Newton- 
OM and Swift-MNOT using the M2 filters. Note 
that these UV data are not corrected for Galactic 
reddening. The XMM-Newton data are adjusted 
to the Swift UVOT data. All fluxes are given in 
units of 10~ 14 W m -2 (=10 -11 ergs s _1 cm~ 2 ) as 
listed in Tabled 




Fig. 3.— Swift XRT and UVOT-UVW2 light 
curves of Mkn335 between 2007 May 17- June 02. 
The upper panel shows the XRT count rate, the 
middle panel the XRT hardness ratio (see § 13.1. 2p , 
and the lower panel displays the UVOT-UVW2 
light curve. 



Fig. 2.— 2007 XMM-Newton EPIC pn light curve 
in the 0.2-10.0 keV band (left) and Swift XRT 0.3- 
10.0 keV light curve of the second Swift o bserva- 
tion on 2007 May 25 (|Grupe et al.ll2007bl ) on the 
right panel. The starting times of the observations 
are 2007 July 10 17:09 UT and 2007 May 25 00:01 
UT, respectively. 
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Fig. 7. — Fits to the 2006 (high-state) and 2007 (low-state) XMM-Newton pn spectra with reflection models. 
The left panels display the fits using model A (i.e. two power law components both with associated reflection) 
and the right panels for model B (i.e. two power law components, of which only one (the flatter) with 
associated reflection). The upper panels display the fits to the high-state, the middle panels the low state, 
and the bottom panels the residuals of the model fits. The line spectra (mekal + emission lines) are in black, 
power laws are shown in red, both reflection components are in blue, and the total model is green. The 
spectral parameters of the joint fits in the 0.3-12.0 keV band are given in Tabled 
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Fig. 9. — RGS spectra of Mrk 335. The upper panel displays the RGS spectra fitted by a single power law 
model. Clearly the residuals reveal several soft X-ray emission lines. The lower panels display the theoretical 
positions of some expected emission lines which are marked in the spectrum at different wavelengths ranges. 
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Fig. 10. — Optical spectrum of Mrk 335. The upper panel displays the whole optical spectrum taken on 
2007 September 07. The lower two panels show the blue and red regions of the highly ionized coronal Fe 
line in the Fell-subtracted optical spectrum of Mrk 335. 
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Fig. 4.— XMM-Newton EPIC pn 0.3-12.0 keV 
(black) and MOS 0.3-10.0 keV (red + green) spec- 
tra of Mrk 335 fitted with a broken power law 
with ax,soft=l-89, orchard = -0.19, £ b r 0a k=l-73 
keV as listed in Table[5] The absorption column 
density is fixed to the Galactic va lue (3.96 x 10 20 
cm -2 ; Dickev fc Lockman 1990l ). 




Energy (keV) 

Fig. 6. — Difference spectrum between the 2006 
high-state and 2007 low-state XMM-Newton spec- 
tra fitted in the 2-5 keV range by a Galactic ab- 
sorber power law with ax=l-3 and extrapolated 
to lower and higher energies. 




Energy [keV] Energy [keVJ 

Fig. 5.— XMM-Newton EPIC pn 2.0-12.0 keV 
spectrum of Mrk 335 with a power law and a broad 
Fe line at 4.9 keV (left panel) and a power law with 
partial covering absorber (right panel). 
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Fig. 8.— Combined XMM-Newton pn and Swift 
BAT spectrum fitted with a power law with partial 
covering absorber. 
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